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Ph! abstract 
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We present ACS/HST coronagraphic observations of HD 100546, a B9.5 star, 
103 pc away from the sun, taken in the F435W, F606W, and F814W bands. 
• Scattered light is detected up to 14" from the star. The observations are consis- 

tent with the presence of an extended flattened nebula with the same inclination 

2 ■ as the inner disk. The well-known "spiral arms" are clearly observed and they 

trail the rotating disk material. Weaker arms never before reported are also seen. 

03 . 

The inter-arm space becomes brighter, but the structures become more neutral 
in color at longer wavelengths, which is not consistent with models that assume 
that they are due to the effects of a warped disk. Along the major disk axis, the 



> 

colors of the scattered-light relative to the star are A(F435W — F606W) ps 0.0- 
0.2 mags and A(F435W — F814W) ~ 0.5-1 mags. To explain these colors, we 
explore the role of asymmetric scattering, reddening, and large minimum sizes on 
ISM-like grains. We conclude each of these hypotheses by itself cannot explain 
the colors. The disk colors are similar to those derived for Kuiper Belt objects, 
suggesting that the same processes responsible for their colors may be at work 
here. We argue that we are observing only the geometrically thick, optically 
thin envelope of the disk, while the optically thick disk responsible for the far-IR 
emission is undetected. The observed spiral arms are then structures on this 
envelope. The colors indicate that the extended nebulosity is not a remnant of 
the infalling envelope but reprocessed disk material. 
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Introduction 



Resolved images of circumstellar disks in optical scattered light provide a unique way 
of establishing the spatial distribution of circumstellar material. They allow for the highest 
spatial resolution of all techniques available and are sensitive to light scattered by small 
particles regardless of their temperature. However, at optical wavelengths one has to contend 
with the brightness of the host star and the use of a coronagraphic mask becomes a necessity. 
Scattered-light observations complement thermal observations, and constrain models based 
on spectroscopic data only. 

In this paper, we present coronagraphic optical images in which the circumstellar mate- 
rial around the B9.5V Herbig Ae/Be star HD 100546 is resolved in scattered light. HP 10054 6 
is a single star, with d=103^g pc and reddening of Ay =0.28 (Ivan den Ancker et al.lll998l ). 
located at the ed ge of the Lynds dark c loud DC 296.2-7.9. Isochrone fitting using the mod- 



els de veloped by iPalla fc Stahlerl (119931 ) suggests its age is >10 Myr (Ivan den Ancker et al. 



19981 ) . The age cannot be much larger than this, as the sy stem presents evidence of strong 
accretion in the H Q lines ( jVieira. Pogodin. &: Franco! Il999l ). Its fractional excess luminosity 
is Lm/L* = 0-51 (Ltr is the excess infrared luminosity and is the stellar luminosity, see 
Meeus et al.ll200lh . 



The circumstellar material was first resolved from the ground in the J and K s bands, us- 
ing the Adaptive Optics Near Infrared System ( ADONIS) mounted on the Europea n South- 
ern Observatory's 3.6 m telescope at La Silla (IPantin. Waelkens. fc Lagagd l2000l ). Those 
coronagraphic observations revealed nebulosity at angular distances ranging from 0."4 to 2" 
from the star. The authors concluded that the observ ations were cons i stent with the presence 
of a disk with an inclination 50° ± 5° from edge-on. lAugereau et al.l (120011 ) performed Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) coronagraphic observations at 
1.6 /im and detected scattered light from ~0."5 to 3. "8 from the star, as well as significant 
azimuthal asymmetries in the circu mstellar disk. The observations had a spatial resolution 
of ~0."2 (IMobasher fc Royd 120041 ) . The highes t spatial reso l ution images to date (~0."1 
per resolution element) have been published by iGrady et al.l (120011 ). These are broadband 



1 For the astro-ph version, the images have been reduced in size and resolution. Please contact D. Ardila 
for a full resolution PDF. 
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coronagraphic observations taken with the Space Telescope Imaging Spectrograph (STIS) 
that probe angular distances as small as ~0."8 from the star (at selected azimuthal angles). 
The observations reveal a complex circumstellar structure which the authors interpret as an 
inclined disk surrounded by an envelope of material extending up to ~10" from the star. In 
addition, the disk presents bright and dark lanes, reminiscent of spiral structure. Thermal 
observations at ~10 /zm using nulling i nterferometry reveal that the disk has a radius of 30 
AU at that wavelength (ILiu et al.l 120031 ). Those observations also suggest the presence of an 
inner hole. 

The spectral energy distribution (SED) of this system has been intensively studied. The 
star is an example of a class of Herbig Ae/Be objects that show substantial silicate emission 
at 10 /im and for which there is strong emission from the mid-infrared (mid-IR) relative t o 
the near- infrared (near-IR). The se have been called group la sources by iMeeus et al.l (120011 ). 
AB Aurigae (4 Myrs, A0V, see van den Ancker et al. 1997 ; DeWarf et al. 2003T) is another 
group la source, often mentioned in the same context as HD 100546. IMeeus et al.l (120011 ) 
suggested that the SED for group la sources could be explained by a geometrically thin, op- 
tically thick disk responsible for the far-infrared (FIR) emission, as well as a flared optically 
thin component responsible for mid-IR emission. However, they did not develop quantita- 
tive models to c onfirm if this configuration could indeed reproduce the SED. More recently, 
Vinkovic et al.l (120061 ) concl uded from the surfa ce brightness profiles derived from the STIS 
coronagraphic observations (IGrady et al.ll200ll ) that the configuration of the material close 
to the star is that of a geometrically t hin, optically thick di sk surrounded by a flat, geomet- 
rically thick, optically thin envelope. IVinkovic et al.l (120061 ) used this model to successfully 
reproduce the near-IR SED (< 6/im). 

For HD 100546, the 3 /i m excess is even small er and the mid-IR excess larger than 
in other Herbig Ae/Be stars (jBouwman et al.l 120031 ). This is significant because, histori- 
cally, the source of the (generally) strong 3 /im "bump" in Herbig Ae/Be stars has been 
difficult to explain, as it is not obvio usly produced by an optically thick disk. For example, 
Hartmann. Kenyon. fc Calvetl ( 19931) have suggested that i t is due to emission by an infalling 
envelope. Recently, Dullemond. Dominik. fc Natta (j2001 ) showed that the puffed- up inner 
edge of a Chiang- Goldr eich disk (a flared disk wi th well-mixed dust and gas and a super- 
heated atmosphere, see IChiang fc Goldreichlll997l ) subjected to direct irradiation from the 
star will emit in the near-IR. If this inner rim is the source of the 3 /im bump, the low 3 /im 
emission from HD 100546 sugge sts a modest rim (outside the bounds of the self-consistent 
model by lDullemond et al.l 120011 ) and a shallow dust surface density f or the inner disk. With 
a model of these characteristics on a 400 AU disk iDominik et al.l (120031 ) were able to fit 
the system's SED, although the fit underpredicts the submillimeter flux by almost an or- 
der of magnitude. Because of the small 3 /im excess, lHarker et al.l (120051 ) were able to fit 
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the SEP with an u nmodified Chiang- Goldreich disk (no puffed-up rim), 150 AU in radius 
Harker et al.l (120051 ) did not include any submillimeter measurements in th eir SEP model, 



which may explain some of the difference in radius with the models from iPominik et al. 
j2003h . 



A completely different model was developed by iBouwman et al.l (120031 ). who fitted the 
SEP by a set of three optically thin spherical shells. The purpose of their model was not to 
determine the geometry of the emission but the mass distribution of the dust as a function 
of temperature. They concluded that the low near-IR emission can be explained by a mostly 
empty region (an inner hole) within 10 AU and speculated that the hole had been created 
by a protoplanet. 

Evidence for a hole close to the star has als o emerged from spectroscopic data. Ultravi- 
olet observations using STIS (IGrady et al.ll2005l ) reveal a deficit of molecular hydrogen close 
to the st ar, which the authors int e rpret as a central gas and dust cavity, 13 AU in radius. In 
addition lAcke fc van den Anckerl ( 2006 ) have used the 6300 A line of [OI] to conclude that 
there is a 4.3 AU-wide gap in the disk, 6.5 AU away from the star. They argue that this is 
consistent with the presence of a 20 Mj object at that position. 

T he mid-IR spectrum of HP 100546 is similar to that of Comet C/1995 01 (Hale- 
Bopp) iMalfait. Bogaert. &: Waelkensl (119981 ) in that it shows strong emission in crystalline 
silicates. In this respect it resembles t hat of the optically thin debris disk around the K0 
dwarf HP 69830 (jBeichman et al.ll2005l ). Models of the material around HP 100546 suggest 
a relatively large fraction of small (0.1 /xm ) forsterite and large (1.5 /xm) enstati te grains 
( iMalfait et al.lll998l ; lvan Boekel et al.ll2005l ). The models by lBouwman et al.l (120031 ) indicate 
that large amorphous silicate grains and small forsterite grains have different temperatures, 
unlike the case in Hale-Bopp, although th ey overlap in rad i al dis tance from the star. The 
mineralogical composition of the dust led IBouwman et al.l (120011 ) to suggest that thermal 
annealing may not have been the origin of the forsterite and that it may result from the 
destruction of a differentiated larger parent body. All these spectroscopic models indicate a 
highly evolved disk both in terms of particle size and crystalline fraction. 

Optical and ultraviolet observations ind icate that, in spit e of the inner hole, the star 
is actively accreting and has a stellar wind (jVieira et al.l 11999c Peleuil et al.l I2004T) . While 



hot m olecular hydrogen has been detected close to the star (ILecavelier des Etangs et al. 



2003 ), the cold ga s component has not been detected (ILoup et al. 



199 



0; 



Nvman et a 



1992 



Wilner et al.ll2003l ). From measurements of the 1.3 mm dust continuum iHenning et al.l (119981 ) 
concluded that there are 240 M ffi of dust within a region 11" x 11" in size, which would 
imply 100 x larger gas mass, if the mixture has ISM prop ortions. This i s cert ainly an upper 
limit as the non-detection of the HCO + line at 89 GHz by I Wilner et al.l (120031 ) suggests that 
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the gas may be strongly depleted. 



The lack of detectable gas, the presence of large particles, and the evidence for an 
inner hole suggest that HD 100546 is an evolved pre-main-sequence star, a "transitional" 
disk perhaps on the way to becomin g a more collisionally ev olved object like HD 141569 
(IClampin et al.l 120031 ) or /3 Pictoris (IGolimowski et al.ll2006l ). In this, it joins a class of 
objects with signs of evolution in their inner disk (flux deficits in the nea r IR relative to th e 
median SED of the class) which have been identified in the last decade (INatta et al.ll2006l ). 
However, the exceptionally high fraction of crystalline silicates and the characteristics of its 
SED suggest that it may be in a class by itself. 

Because of its relevance in this poorly understood tra nsitional phase, we targeted it for 
observation using the Advanced Camera for Surveys (ACS; iFord et al.ll2003l . |Pavlovskill2006l ) 
on board the Hubble Space Telescope (HST). In this paper, we present F435W, F606W, 
and F814W coronagraphic images of HD 100546's circumstellar environment. These images 
reveal the environment between 160 and 1300 AU from the star with unprecedented spatial 
resolution, suppression of the wings of the point spread function, and photometric precision. 
In addition, they allow us to measure for the first time the color of the circumstellar dust. 
The scattered-light observations presented here sample the small dust population, far (>160 
AU) from the star. They are then complementary to spectroscopic observations and SED 
modeling that focus primarily on the inner disk. 



2. Observations and Processing 

In this section we present a summary of the reduction procedure, which involves the 
subtraction of a point spread function (PSF) from the target star. The reduction follows the 
standard procedures common to other coronagraphic observations. The resulting images were 
then deconvolved using a Lucy-Richardson algorithm. Those deconvolved images form the 
core data used for the analysis. The science- quality images on which most of the conclusions 
in this paper are based are available from the Journal's Web site. 

These observations were conducted as part of the guaranteed observing time awarded 
to the ACS Investigation Definition Team (proposals 9987 and 9295). They were performed 
using the coronagraphic mode of the ACS High- Resolution Channel (HRC), with t he 0" .9 oc 



cult in g spot, which suppresses the stellar PSF wings by factors of the order of ten (jPavlovski 



20061 ). It is possible to improve the contrast by subtracting a star of the same spectral 
type as the target, which suppresses the remaining PSF wings. To do this subtraction, we 
observed HD 129433 as a PSF reference star. The star was chosen because it is bright, has 
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similar color (V=5.73, B-V=-0.01. ICrawfordlll965t ). 

The observing log is shown on table [TJ The images were taken with the F435W (Johnson 
B), F606W (broad V) and the F814W (broad I) filters. The HRC has a pixel scale of 25 mas 
pixel -1 , and a coronagraphic field PSF FWHM of 50 mas in the F435W passband, 63 mas in 
the F606W passband, and 72 mas in then F814W passband. For each band and each star, a 
short direct exposure was followed by a coronagraphic exposure. In each band, HD 100546 
was observed at two different telescope roll angles. The process of PSF subtraction produces 
image artifacts and the two orientations help distinguish between these and real features. 



2.1. Reduction 

The initial stages of the image reduction (i.e., subtraction of bias and dark frames, and 
division by a non- coronagraphic flat field) were performed by the ACS ima ge calibration 



pipeline at the Space Telescope Science Institute (STScI) (jPavlovskil 120061 ) . We started 
our own image reduction procedures on the separate flat-fielded images before combination, 
cosmic-ray rejection, and correction of the HRC's geometric distortion by the STScI data 
pipeline. Coronagraphic images with the same exposure time were averaged and cosmic rays 
were rejected using conventional IRAF toola^l To account for vignetting around the occulting 
spot, the coronagraphic images were manually divided by the spot flats available from the 
ACS reference files page at STScI's Web site. Each spot flat was shifted to the appropriate 
position corresponding to the time of observation. This step is done automatically since the 
CALACS pipeline version 4.5.0. 

The coronagraphic images obtained in each band and roll position were then sky- 
subtracted. The sky was measured by taking the median in the number of counts from 
circles 20 pixels in radius at the four corners of the images. The effect of this step is minimal 
in the overall photometry. After sky subtraction, the coronagraphic images were divided 
by the exposure time. In the longest exposures, light diffracted by the coronagraphic mask 
saturates some pixels in the region within 2" from the star and we replaced those saturated 
pixels by non-saturated pixels from shorter exposures. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.2. A Note on HST's Photometric systems 



The F435W filter is a good approximation of Johnson B, but the other two filters used 
in these observations do not correspond to standard filters. To calibrate our observations and 
compare them with published results, we make use of three different photometric systems: 
an "instrumental" vegamag syst em, a standard vegamag system, and the standard Johnson- 
Cousins system as described by lLandoltl (119921 ) . 



The instrumental vegamag system is described by lSirianni et al.l (120051 ) . In this system 
Vega has zero magnitudes and zero colors in all bands. The transformation from e~/sec to 
magnitudes involves only a zero-point and no color term. We will refer to magnitudes in this 
system with the filter name. For example, the F435W magnitude, etc. This is the system 
most used in this paper. 

If the spectrum of the source is known, it is possible to obtain magnitudes in the 
standard vegamag system. These magnitudes are given in the Johnson B, Johnson V, and 
Cousins I bandpasses, but the calibration is still tied to Vega. Transformations between the 
instrumental and the standard vegamag system are possible using the synthetic photometry 
package SYNPHOT, which is distributed by STScI. We will refer to the magnitudes in this 
system with the subscript Vega. For example, Bv ega , Vv ega , etc. 

Ground based observations usually present measurements in a standard photometric 
system, involving standard filters and a grid of primary calibrators. ISirianni et al.l (120051 ) 
provide color transfo rmations betwe e n the instrume n tal ve gamag system and the standard 
system described by lLandoltl (I1992I ). ISirianni et al.l (120051 ) argue that the transformations 
are strongly dependent on the details of the source spectra and so we use them sparingly. 



2.3. Stellar Photometry 

A measurement of the stellar flux is required for two aspects of the analysis. First, in 
order to subtract the PSF reference star from the target, the former needs to be scaled by 
a multiplicative factor. This factor is the ratio of the brightness of the reference star and 
that of HD 100546 in each band. Very accurate photometry is not crucial for this procedure 
because visual inspection of the subtracted images is enough to determine the scaling factor 
within ±3%. Second, in order to determine the intrinsic color of the circumstellar material, 
we divide the images by the stellar flux. For this, we need to know the stellar photometry as 
accurately as possible. As one of the main conclusions of this paper is that the circumstellar 
material is intrinsically red, we now describe in detail the procedure followed to obtain this 
photometry. 
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For F435W we used the direct (non-coronagraphic) images of HD 100546 and HD 129433 
to measure the number of e~/sec within 5" . While the cores of the stellar PSFs are saturated, 
observations performed with the pa rameter GAIN set to 4 (as are these) preserve the number 
of electrons over the whole image (jGilrilandl 120041 ) . In other words, as long as the aperture 
is large enough to encompass all the saturated pixels, the integrated number of counts can 
be recovered. 

Direct non-coronagraphic observations in the other bands were executed with GAIN=2 
and so the number of counts in saturated images is not conserved. We used the number of 
counts measured in F435W, a model spectrum of the target and the PSF reference star, and 
the measured colors of HD 100546, to predict the number of counts in F606W and F814W. 



For HD 100546, Ide Winter et all fl200lh measured B-V= 0.012 ± 0.005 mags and V- 



van 



den Ancker et al.l Jl998h derived A y =0.28 



mags, as- 



1=0.042 ± 0.01 mags, from which 
suming that the intrinsic color was B-V=-0.075 mags. In order to compare with our data, 
we assumed that the colors in the standard vegamag system are the same as those in the 
standard Johnson- Cousins system. This is a very good assumption for HD 100546, whose 
spectral t ype is B9.5 and so the color- dependent terms in the transformation equations 
are small (jSirianni et al.l 120051 ). As the model spectrum for the target we used the B9V 
star HD 189689, from the Bruzual-Persson-Gunn-Stryker (BPGS) Atlas provided by STScI. 
This star has By ega -Vv Pga =-0.11. We r eddened the spectrum using the standard extinction 
law f Cardelli. Clayton, fc Mathislll989l ) until Bv eg a-Vve g a=0.012 mags, the average measured 
color for HD 100546. The resulting extinction is Ay=0.36 mags. With this spectral model 
and the number of counts measured in the F435W band, we used the calcphot routine within 
SYNPHOT to derive a Bv ega magnitude. The Bv ega magnit ude derived in thi s way_ (table [2]) 



is 0.03 mags less than the average B magnitude derived by Ide Winter et al.l (120011 ). consis- 
tent the difference between the standard system, in which Vega has B=0.03 mags, and the 
vegamag system, for which Vega has Bv eg a = 0.00 mags. From the stellar colors of HD 100546 
given by Ide Winter et al.l (120011 ) we derived Vv ega and Iv eg a- With these magnitudes, the cal- 
cphot routine, and the same spectral model, we obtained the number of counts in F606W and 
F814W. The result of this procedure is that the number of counts in each band corresponds 
to the measured stellar color. 

The BPGS Atlas contains four more dwarf stars with spectral types between B9 and 
AO. For each of them, a different value of Ay is required to reproduce the observed color for 
HD 100546. They range from Ay=-0.077 mags in the case of 58 Aql to Ay=0.18 mags for 
9 Vir (both classified as A0V in the BPGS Atlas). However, following the same procedure 
described above results in changes of the number of counts for each band of <1% with respect 
to those reported in table El The procedure is robust to the assumed spectral model because 
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each model is forced to have to observed colors of HD 100546, independently of its intrinsic 
colors. 

For HD 129433, the number of counts in F606W and F814W was obtained from the num- 
ber of counts in F435W by using as spectral model the unreddened spectrum of HD 189689. 
The PSF reference star has a larger Bv e ga-Vve g a (by 0.1 mags) than HD 189689, and so this 
procedure will systematically underestimate its flux by fh 10%. Accurate photometry of the 
reference star is not crucial for the conclusions of this paper. 



Sirianni et al.l (120051 ) indicates that the best aperture ph otometry using the HR C will 



have 3% errors. For HD 100546, the color measurements from lde Winter et al.l (120011 ) intro- 
duce an additional 2% error (3% error) in the determination on the F606W (F814W) flux. 
For the PSF reference star, we assume that the errors are those of the aperture photometry 
only. These are the errors quoted in table [2J 



2.4. Subtraction of the Coronagraphic PSF 

For each band, the normalized, cleaned, and combined coronagraphic images of the 
PSF reference star HD 129433 were subtracted from each roll of the combined images of 
HD 100546. The positions of the occulted stars were determined within one pixel using the 
central peaks of the coronagraphic PSFs t hat result from th e re-imaging of the incompletely 



occulted, spherically aberrated starlight (IKrist et al.l 120031 ) . Further alignment was done 



visually, by shifting the PSF reference images until subtraction residuals were minimized. 
The uncertainty in the registration of the images along each axis is 0.125 pixels. For the 
normalization constants we performed the following visual, iterative adjustment. Starting 
with the flux ratios implied by the values from table H (0.391, 0.443, 0.474 for F435W, 
F606W, and F814W, respectively), we adjusted the scaling until the subtraction residuals 
were minimized. This occurred for normalization constants 0.383, 0.382, and 0.387 in F435W, 
F606W, and F814W. Visually, the uncertainty of these normalization values is 3% and 
in what follows we will propagate this error linearly (as a systematic error) to estimate 
uncertainties in calculated quantities (Section [33]). The values of the normalization constants 
are consistently lower than implied by the values quoted in table [2], as expected by the 
mismatch between the colors of the PSF reference star HD 129433 and those of its spectral 
model. 

After subtracting the coronagraphic images of HD 129433 from each image of HD 100546, 
we corrected for the geometric distortion in the HRC image plane, using the co efficients of 



the biquartic-polynomial distortion map provided by STScI ( IMeurer et al.ll2002l ) and cubic 
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convolution interpolation to conserve the imaged flux. By using standard IRAF routines, we 
aligned the two rolls in each band, using other stars in the field for reference, and averaged 
the images from each roll. The resulting surface brightness maps (divided by the stellar 
fluxes) are shown in figure [TJ These maps are available from the Journal Web page. 

The mid and bottom rows of figure [1] highlight different regions of the circumstellar 
environment. Differences in the shape of the PSF between HD 100546 and the reference star 
HD 129433 (due to differences in color, telescope focus, or position behind the mask) result 
in subtraction residuals. The radial streaks in the top two rows in figure [JJ are produced 
by differences in the way light is scattered by telescope optics for the two stars while the 
concentric dark and bright ri ngs are due to differences in the diffraction pattern produced 



by the coronagraphic mask (IKristl 1200 ll ). The rings are more noticeable along the NE- 
SW direction, which coincides with the position of the scattering "strip" . This well-known 
artifact of the ACS coronagrap h is a bar-like regio n of increased number of counts, centered 



on the star and ~6" long (e.g. iKrist et al.ll2005bl ). 



The presence of these residuals makes the analysis within m 1.6" of the star impossible. 
Beyond this limit we can use the images in the two separate telescope rolls to distinguish 
subtraction residuals from features associated with the material around HD 100546, as the 
residuals are nearly fixed on the detector's axes. To trace the contours shown in the bottom 
row of figure [1] we perform a 10-pixel median azimuthal smoothing on the middle-panel 
images, followed by a 5-pixel boxcar average smoothing. 

In figure [TJ the difference in the spatial resolution among the three bands is seen as a 
smoothing of the images with increasing wavelength. In the top row, the imprints of the 
larger coronagraphic spot and the occulting bar are visible towards the SW direction (one 
for each roll angle). 



2.5. Red Halo and Image Deconvolution 



For HRC observations in the F814W band, one has to contend with the effect of the 



red halo (jSirianni et al.ll2005l ). This may appear as a diffuse, smooth halo of light, slightly 
asymmetric but fixed with respect to the detectors axes. The halo is composed of photons 
from the core of the stellar PSF which have been scattered to large angles by the CCD 
substrate. The hal o asymmetry can be seen in the direct (non-coronagraphic) images from 
Krist et al.1 (l2005bl ) (See their figure 2). In this section we show that deconvolution of the 
stellar PSF, which partially corrects this effect, changes modestly the slope of the surface 
brightness profiles. 



- 11 - 



To deconvolve the images, we use the P SFs generated by the Tiny Tim software package, 



distributed by STScI (IKrist h Hookll2004l ). The simulated "off-spot" coronagraphic PSFs 
can be generated for an array of spectral types. These PSFs include a simplified model of 
the red halo but do not include the red halo asymmetries mentioned before. We do not 
believe those asymmetries are important in our images: upon subtraction of the two rolls 
no residual asymm etric structures are revealed. Unlike the non- coronagraphic images from 



Krist et al.l (j2005bl ). the coronagraphic images presented here do not produce a very strong 



core to be scattered at large angular distances. 

The deconvolution is performed on the image after subtraction of the reference PSF. 
Therefore, only the occulted star subtraction residuals contribute to the deconvolution pro- 
cess (and then basically as noise) as the effects of vignetting near the occulter are removed 
with the division by the coronagraphic spot flat field. The only effect of the coronagraphic 
setup on emission outside of the focal plane mask is to alter the PSF, which is then defined 



by th e Lyot stop rather than the entrance aperture of the telescope. Simulations (IWhite 



20071 ) confirm that beyond 1.5" from the center of the coronagraphic mask, the off-spot PSF 



does not depend on distance. 

We generated distorted off-spot PSFs in each band, using BPGS model of the AO star 
9 Vir as our model. This is the closest spectrum to our target provided by Tiny Tim . 



We applied the Lucy-Richardson deconvolution algorithm (jLucyl Il974l ; iRichardsonl [1972|) 
as implemented in IRAF to each (single-roll, distorted, PSF-subtracted) image outside a 
circular region 1" in radius, centered on HD 100546. We terminated the computations at 50 
iterations. By this point, we could not detect changes in the FWHM of the field stars, nor 
in the surface brightness profiles of the image. After deconvolution, we corrected for camera 
distortion, rotated the two rolls in each band to a common orientation and averaged them. 
In this way, we obtained a deconvolved image in each band. A comparison of the "before" 
and "after" surface brightness profiles in each band is shown in figure [2J The resulting images 
are shown in figure [31 

Figure [2] shows that the deconvolution has its strongest effect in the F814W images, 
where it results in a sharpening of features within 2." 5 and a modest reduction of the large 
scale surface brightness from 2." 5 to 10" away of HD 100546. As mentioned before, the 
deconvolution process is not valid within l."5 from the coronagraphic mask. In the case of the 
direction shown in figure [2] (corresponding to the SE semi- major disk axis), the deconvolved 
profiles are ~15%, ~5%, and ~3% fainter in average than the non- deconvolved ones, for 
F814W, F606W, and F435W, respectively. These values are similar to those along other 
directions. Therefore, the A(FA35W — F814W) color of the deconvolved images (defined in 
section |3T4"|) is ~0.12 mags less (i.e. not as red) than the color of the non- deconvolved ones. 
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As we will show in section [331 this is comparable to the errors in the color. So in terms of the 
large scale behavior of the color it does not make much difference if the analysis is carried 
out in the deconvolved or in the non-deconvolved images. 

Figure [3] shows the deconvolved equivalent of figure [TJ After deconvolution, the features 
in the circumstellar environment become sharper, and the images become noisier. The 
amplified, correlated noise in the deconvolved images is characteristic of the Lucy-Richardson 
algorithm and it is a consequence of the algorithm's requirements of a low-background signal, 
positive pixel values and flux conservation in local and global scales. The higher noise levels 
require that the images be further smoothed before calculating the contours. To measure 
them, we smooth the images as in figure [1] and then convolve them with a gaussian function 
having FWHM=10 pixels, which acts as a low-pass filter. It is in these heavily smoothed 
figures that the contours are traced. 

The deconvolved F435W, F606W, and F814W surface brightness maps are available 
from the Journal Web page. 



2.6. Color-combined images 

In order to grasp the overall appearance of the system, we created color composites of 
the images by mapping the bands in three color channels: F435W in Blue, F606W in Green 
and F814W in Red. Before combining them, each image was divided by the stellar flux in its 
band. To produce the color composites, the images were scaled by the inverse hyperbolic sine 
of the average in al l three bands, using the softening parameter (3 = 0.005, as described in 



Lupton et al.l (12004 ) . This scaling makes the apparent color independent of the brightness. 



In effect, the scaling is linear for very faint regions and logarithmic for brighter ones. 

The results are shown in figure HI Regions of extreme colors, like the circular purple 
features close to the coronagraphic mask, indicate places in the image for which data from 
at least one of the bands are missing. These occur because the position of the coronagraphic 
subtraction residuals is wavelength dependent. The deconvolution process (figure HI right) 
makes the radial streaks resulting from PSF mismatch very evident. 



3. Results 



Figures dj [31 El show the well known circumstellar structure of HD 100546 (IGrady et al. 

20011 ) but with an unprecedented degree of clarity. Within l."6 the images are dominated by 
subtraction residuals. From l."6 to 10" from the star, the total scattered-light relative fluxes, 
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measured in the deconvolved images, are (Lcircum/L*)x = 1.3 ±0.2 x 10 3 , 1.5 ±0.2 x 10 3 , 
and 2.2 ± 0.2 x 10~ 3 , in the F435W, F606W, and F814W bands respectively. These values 
are ~ 200 — 400 times less than the total IR excess. 

In this section we describe the main observational results associated with the general 
characteristics of the scattered light distribution around HD 100546. 



3.1. Photometry of Field Stars 



HD 100546 has a galactic latitude of -8.32 degrees, which suggests that the other point 
sources in the images are unassociated field stars. In table [3] we list the standard photometry 
for the brightest stars within ~6" of HD 100546. See also figure [51 We only list those 
objects that are observed in both rolls (and are outside the occulting bar in both) in at least 
one band. To compare with known stellar colors, the ph otometry is given in the standard 
Johnson- Cousins system, using the color corrections from ISirianni et al.l (120051 ). 



The largest V-I color is ~ 2. If the object is unreddened, this color would suggest that 
it is an ea rly M dwarf. Such a star would have an I magnitude of ~ 12, much brighter than 
observed (lLeggettlll992l ). All other objects have bluer colors, indicating earlier spectral types 
and brighter predicted apparent magnitudes. Therefore, all objects in table [3] are farther 
away than HD 100546. 



3.2. General Morphology 



Overall, the images show a broad band of light wrapping around the SW side of the 
star (figu re HI top), reminisc ent of the configuration of the scattered light images around AB 
Aurigae (IGrady et al.l 119991 ). A red glow can be seen on the NE side of the target both in 
the deconvolved and non-deconv olved images, similar t o the strong F160W scattered light 
signal found closer to the star by lAugereau et al.l (120011 ). 



Interpreting the inner circumstellar surface brightness as due to an inclined disk, we 
use the deconvolved images to fit elliptical isophotes to the images, with semi-major axes 
ranging from l."6 to ~ 2". At distances larger than 2" the contours become "boxy" because 
of the bright SW feature (structure 2, see section 13.31) and the scattered light cannot be fit 
with an ellipse. For all the bands, the position angle (PA) of the semi-major axis is ~ 145 
degrees. From the ellipticity of the isophotes we derive an inclination angle of 42 degrees 
from face-on. There are no significative differences in these values among the three bands. 
This is an observational determination of the inclination and PA based on isophotal contours, 
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and ignores the effect of anisotropic scattering of the stellar light by the dust grains. These 
results are summarized in table HI 

Table H] also compares the disk geometry derived in this work with other published 
results. All disk inclinations are consistent, but t he position angles of the major axes ar e 
only marginally so. Our res ults differ from th ose bvlPantin et al.l (120001 ). iGrady et al.l (]200ll ). 



and lAugereau et al.l (120011 ) at t he 2a level. iPantin et al.l (120001 ) derived the PA by using 
symmetrized imag e of the disk. iGradv et all (bpOlh do not describe the procedure used to 



determine the PA. lAugereau et al.l (1200 ll ) determined the PA by performing isophot fitting 
within 3" . These differences may reflect the fact that a measurement of the ratio of the axes 
of an ellipse is less sensitive to uncertainties than measurements of a single axis. Furthermore, 
there may be differences between the value of the scattering phase function in the optical 
and the near-IR. 

For the F814W images, the contours within 2" (bottom row of figures [1] and [3]) are closer 
together at the SW side than at the NE side: the image brightness changes faster in the 
former than in the later. The contours in the other bands are not as well defined at these 
distances. 

In the sections that follow and for the purposes of the morphological description, we 
divide the circumstellar environment in three regions, where the distances are measured 
along the disk major axis: the inner disk (from l."6 to 3"), the mid-disk (from 3" to 8") 
and the extended envelope (beyond 8"). Each region has its own dynamical and/or radiative 
characteristics. The inner disk and the material under the coronagraphic mask, comprise 
the region believed to be responsible for most of the far-IR excess. 



3.3. Structures in the Disk 



The well-known "spiral arms" of HD 100546 (e.g. IGrady et al.l 120011 ) are clearly seen 
beyond l."5 to ps3" from the star, as elongated structures at the NW, SW, and SE (figure E]). 
The simpler mask shape of the ACS/HRC coronagraph compared to the STIS one, allows 
for a more complete description of the environment. In the rest of the paper, we refer to 
these features as structures or spiral arms, interchangeably. 

Structure 1 is a narrow band of brightness visible between 2." 3 and 2." 8 from the star, 
and figure [6] suggests that it may continue on the NE side of the disk. Structure 2 describes 
a broad arc at 2", and is separated from the inner part of the emission by a dark region 
~0."5 wide. However, the exact extent of the dark lane is hard to judge because of the 
subtraction residuals along the minor axis. Structures 1 and 2 are separated by a region 
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that is half as bright as the maximum of structure 2. Structure 3 extends from l."5 to 
2." 5 from the star, and it seems more open than structure 1. Because of this, structure 1 
seems to wrap itself around the NE side of the disk, while structure 3 stabs the SW side. 
Structures 1 and 3 are superficially reminiscent of galactic spiral arms. The morphology of 
structure 2 is reminiscent of that produced when starlight is scattered fro m the disk back 



side t owards the observer (see for example, the observat ions of GM Aur in ISchneider et al. 



20031 ). However, both in observations and simulations (jWhitney fc Hartmannl 119921 ). such 
structure is generally narrower and it wraps around the disk over a larger angle than seen 
here. 



Observations of the [01] 6300 A line by lAcke fc van den Anckerl ( 20061 ). reveal blueshifted 



emission from the SE of the disk and redshifted emission from the NW part, suggesting that 
the inner part of the disk is rotating counter-clockwise. Their observations trace gas at 
distances <100 AU from the star. Assuming that the disk rotates in the same direction at 
larger distances from the star and that the SW side is oriented towards the Earth (section 
I3.4.2p . those spectral observations indicate that the structures are trailing the direction of 
rotation. 

Before deconvolution, the behavior of the structures with wavelength is tangled with 
the changing resolution of the telescope at different passbands (figured]). The deconvolved 
images (figure [3]) reveal no obvious morphological differences in each band, for each of the 
structures. The contour levels traced on the deconvolved images show that the space between 
the inner part of the disk and the structures becomes brighter at longer wavelengths, while 
the value of the peak brightness does not change much between bands. This color behavior 
is discussed in the section 13.41 

The circumstellar disk is illuminated by the starlight, which decays as the inverse dis- 
tance squared from the star. An appropriate correction for this effect is only possible for 
face-on disks. We therefore deproject the disk by the inclination, divide by the Henyey- 
Greenstein scattering phase function with g = 0.15 (derived in section 13.4.21) . and multiply 
every pixel by the square of the projected distance to the star. This procedure will reveal 
the correct geometry and brightness of the circumstellar material only in the case of a disk 
in which every dust particle is illuminated by the full stellar radiation field. Such is not the 
case for HD 100546, as every reasonable model of the SED indicates that some part of the 
disk is optically thick. However, independently of the optical depth, the procedure helps 
to reveal weak disk features and clarifies the disk structure, although the brightness of the 
structures and their distance to the star, particularly along the minor axis, will be incorrect. 

The result is shown in the top left panel of figure [7] for the non-deconvolved images in 
the F435W band (images obtained from other bands and from the deconvolved data look 
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similar). Because the hemispheric subtraction residuals are circular in the original image, 
they appear elliptical in the deprojected one. In order to accentuate the sharp structures, we 
subtract from the deprojected image a smoothed version of itself (figure [7J top right). The 
smoothed image has been convolved with a gaussian function having a 30-pixel FWHM. The 
subtraction of this smoothed image from the deprojected one leaves only the sharp features. 
An unsharp masking te chnique similar to this has recently been used to clarify the arm 



structure of our galaxy (jLevine. Blitz, fc Heilesl 120061 ) . 



This unsharp masking also sharpens the subtraction residuals present in the images, 
and so identification of real features requires care. In the bottom right panel of figure [7J we 
have marked those features that are observed in all three bands. The original structures 1, 
2, and 3 are traced with solid lines, and marked as la, 2a, and 3a, while new structures are 
traced with thick dashed lines. 

The brightest new features are structure lc, an inner arm at ~ 2" that seems to emanate 
from la, and structure Id, a feature behind the disk. The deprojected images show that 
structure 2a curves toward the star, like another spiral arm. The la+lb arm is ~5" long, 
while 2a+2b and 3a+3b are ~4" long. On the deprojected, unsharp masked image, arms 
la and 3a have the same opening angle with respect to the disk and appear like symmetric 
images of each other. Together with arms 2a and Id, the overall configuration is that of a 
four-armed spiral. The unsharp masking traces only the sharp ridges, and so the features 
may extend beyond these limits. 

Structures analogous to these have al so been detected aroun d AB Aur, based on ground- 



based coronagraphic images in the H-band lFukagawa et al.1 (120041 ). although they do not seem 



as well ordered. In those observations, at least one of the arms appears split, as structure 
1 is split between la and lc. The arms around AB Aur are also believed to be trailing the 
disk. 



3.4. Surface Brightness Profiles 

Figures [8] and [9] show the relative surface brightness profiles (SBPs) and colors, along the 
major and minor axes of the disk. The SBPs are obtained from the deconvolved images by 
taking median values along rectangular strips (0."25 wide) centered on the star for each band. 
In addition, each profile has been median-smoothed with a boxcar function 10 pixels long. 
The bottom plots show the colors measured with respect to the star: A(FA35W — F606W) = 
-2.5\og(SBP F435W /SBP Fm6W ), A(F435iy - F8UW) = -2.5 log(SBP F43SiW /SBP F&uw ). 
An aid to the interpretation of the profiles is provided in figure [TDJ 
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3.4-1- Sources of Error and Comparison with Previous Measurements 

There are three main sources of error in the SBPs: photometric errors due to the 
uncertainty in the stellar flux (which varies between 3 and 5%, as given in tableland result in 
errors small compared to those from other sources), photometric errors due to the uncertainty 
in the flux ratio between HD 100546 and HD 129433, and errors (variability) within the 0."25 
segments. In figures [8] and [9] the total la errors are shown for some representative points. 
For these profiles, the 3cr detection limits are 2 x 10~ 7 /arcsec 2 , corresponding to Vv eg a=23.4 
magnitudes/arcsec 2 . 

The uncertainty in the flux ratio produces errors which scale with the brightness of the 
circumstellar material. We propagate them linearly, as they are systematic and not random. 
For the major axis, they range from 10% for distances <8" to 5% farther out. For the minor 
axis they range from 20 to 40%. The errors are larger along the minor axis because that 
axis coincides with the position of the scattering strip (see section I2T41) . Finally, errors in the 
median within the 0."25 segments are a measure of the variability within the strips used for 
the calculation of the SBPs, and range from ~5%, 2" from the star, to ~10% at 12" for the 
major axis or ^40% at 8" for the minor one. 

In summary, errors in the photometry (introduced by errors in the normalization con- 
stant) dominate the error budget at small distances and uncertainties due to variability 
within the median strip dominate at large distances, with a turnover point at 6-8" from the 
star. 

In addition, there are regions of the profiles strongly affected by the systematic mis- 
matches between HD 100546 and the PSF reference star HD 129433. For example, the "dip" 
in the SE semi-major axis in F606W at l."5 is due to over-subtraction of the PSF at that po- 
sition, as it is obvious from figure [3j Within « l."6 the SBPs in all directions are dominated 
by subtraction residuals. 

As mentioned in the introduction, previous coronagraphic obse rvations of th e circu mstel- 



lar environment are available. Over the unfiltered STIS passband, iGrady et all (120011 ) mea- 
sured a total count rate of 126 cts/sec at 2", on 9 pixels, which implies 5382 e~/sec/arcsec 2 . 
With this measurement, the reddened spectrum of HD 189689 from section [2~3| and calcphot, 
we predict that the count rate in the F606W passband should be 5 x 10 3 e~/sec/arcsec 2 , or 
1 x 10~ 4 per arcsec 2 (after dividing by the stellar flux in the F606W passband). This point is 
shown in figure [HJ We repeat the analysis for a point at 8" . Both points are broadly consis- 
tent with our observation s . To c ompare with the NICMOS coronagraphic images we use the 



fact that lAugereau et all (120011 ) measured 16 mag/arcsec 2 at 2." 5 from the star, which im- 



plies 1.7x10 4 /arcsec 2 , at the F160W band. To obtain this number we used the magnitude 
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zero-points reported by lAugereau et aJj (120011 ) and the reddened spectrum of HD 189689 as 
a spectral model, which implies that the F160W flux density of HD 1 0546 is 2.5 Jy. We 
also plot the NICMOS measurement at 3." 5 for which lAugereau et al.1 (120011 ) measured 18 
mag/arcsec 2 . The NICMOS measurements imply A(F435W - F160W) =1.5 at 2. "5 and 
A(FA35W — F160W) =1.4 at 3." 5. As before, uncertainties in the stellar model trans- 
late into uncertainties of ~ 4% in t he calculated fl ux. C omparison with the ground-based 
ADONIS measurements reported by iPantin et al.1 (120001 ) yields results that are difficult to 
understand. They measured 0.3 Jy/arcsec 2 at 0."5 from the star in the J-band. This is 
one order of magnitude larger th an the F160W NICM OS measurement at that point 11 
mag/arcsec 2 or 0.04 Jy/arcsec 2 ). ISchutz et al.1 (120041 ) has shown that variations in the PSF 
shape due to atmospheric fluctuations may result in false disk detections, when using ADO- 
NIS paired with co ronography. On the other hand, the PA and inclination reported by 



Pantin et al.1 (120001 ) are roughly consistent with those derived using other instruments. 



3.4-2. Morphology of the Nebulosity 

The top panels of figure [8] show that at 8" the F814W SBP drops by a factor of ~2. The 
drop is very abrupt along the SE direction. Beyond this distance, the A(FA35W — F814W) 
and A(FA35W — FQ0QW) colors become indistinguishable from each other. Consistent with 
the measured inclination, there is a decrease along the SW minor axis at 6" -7". This is 
why we define the mid-disk as the disk between 3" and 8" from the star, and the extended 
envelope as material beyond 8". 

Figure [TU] shows that the nebulosity is slightly more extended to the NW side than to 
the SE. Along the NW side, the major-axis SBPs show scattered light to the edge of the 
detector, ~ 14". Along the SE side the emission drops at 12", from 6 x 10 -7 /arcsec 2 to below 
the detection limit. Again, this is consistent with the decrease along the NE semi- minor axis 
at ~ 9". In other words, the extended envelope is a flattened nebular structure ~12" - 14" 
in radius, seen with an inclination similar to the inner disk. 

In figure [TT] we compare the profiles along each semi-axis to each other. The semi-major 
axis profiles along each direction are very similar to each other, while the profiles along the 
semi-minor axis are very different. Along the major axis, the main differences between the 
two sides are in the presence of structure 1 at 2." 8, structure 3 at 2." 3, and field star ^ 8 at 
5.5". Along the minor axis, the amount of scattered light in all bands is larger in the SW 
side than in the NE side, until the sharp decrease of the SW side at 6" -7". The ratio of 
the scattered light in the SW side to that in the NE side, in each band, between 2" and 7" , 
presents strong local variations, with values oscillating between 1 and 4. The ratio averaged 



-19- 



over all bands in this range is 2.3. For larger distances, the ratio at all bands is consistent 
with one. 

If the brightness asymmetry between the SW and NE sides is interpreted as being the 
result of forward scattering, it implies that the SW is oriented towards the Earth. Assuming 
that the scattering phase fu nction can be parameterized as a Henyey-Greenstein function 
( iHenyey fc Greensteinlll94ll ) the average value of 2.3 for the SBPs ratios implies a scattering 
asymmetry factor g fa 0.15. A rati o of 4 would imply q pa 0.23. For comparison, ISM grains 
have g~ 0.4 to 0.6 in the optical (jWeingartner fc Draindl200ll ). In other words, the grains 
around HD 100546 scatter light more isotropically than grains in the ISM. This determination 
of g is not strictly correct as the procedure is only valid for optically thin dust distributions. 



3.4.3. Power-law Fits 



In any direction, the SBPs are considerably more complex than any previously published. 
For the purposes of comparison with previous works, we derive the power-law fits (see table 
[5]) for the SE Major axis. The profile can be broken in four parts: l."6 to 5", 5" to 8", and 
8" to 12". To do the fit, we interpolated over the signature of structure 3. The results of 
the fit for F435W are shown in figure (TTJ 

While it seems clear that the profile to 5" can be described by a power-law, such de- 
scription is less warranted beyond this limit. In particular, there seems to be a transition 
region from 5" to 8" , before a flatter edge at longer distances. Within the errors, the fits are 
the same across wavelengths. 



Power-law fits to SBPs were derived by lGrady et al.l (120011 ) based on STIS observations, 
using a different PA for the disk, different width for the evaluation of the SBP, and different 
distance ranges fo r the power-laws. In table [5] we also derive the power laws using the 



Grady et al.l (120011 ) ranges and tracks, but based on our data. For the exponents, they derive 
-3.1 ±0.1 between 1" and 2." 7, -3.1 ±0.1 between 2." 7 and 5", and -2.2 ±0.2 between 5" 
and 8" . Within the errors, their results are consistent with ours except in the range between 
2." 7 and 5", were our slopes are significantly steeper. The origin of this discrepancy cannot 
be ascribed to the differing observation bands, as the power- law indices derived in all our 
filters are alike within the errors. From NICMOS F160W data. lAugereau et al.l (120011 ) derive 
-2.92 ± 0.04 between 0."5 and 2."5 and -5.5±0.2 from 2."7 to 3."7. Between 2."7 to 5."0, 
our measured index is between those of STIS and NICMOS. The difference with NICMOS 
may be due to the differing optical depths of the material at different wavelengths. 



Beyond 5", 



Grady et al 



(120011 ) fitted a power-law index of -2.2 along the major axis 



-20- 



of the disk. Based on this fit, and the models from IWhitney &: Hartmannl (119931 ). they 
concluded that they had detected an in falling envelope in the star. The applicability of the 
models by lWhitney fc Hartmannl (119931 ) presumes that the material is optically thick beyond 
500 AU from the star. By using the same position angle and stri p size to mea s ure th e major 
axis, we can roughly reproduce the power-law dependence from iGrady et al.l (120011 ). within 
the errors (See table [5]). However, at large distances from the star, the material will be 
optically thin along the line of sight. The index of the power-law then reflects the surface 
density of scatterers, and not the shape of the sca ttering surface. As such, the use of the 
models developed by IWhitney fc Hartmannl (119931 ) to explain the large scale power-law is 
unwarranted. 



3.5. Disk Color 

The SBPs show that along the principal axes of the disk the circumstellar material is 
brightest in F814W than in the other bands. Furthermore along the major axis, A(F435W — 
F814W) increases from 0.3 at 2", to ~ 1 at 5", indicating that the dust is between w 1.3 
and 2.5 times brighter in the F814W band than in F435W. After 8" there is a sharp decrease 
in the value of A(FA35W — F814W). The noise precludes an accurate determination of 
the value, but at least from 8" to 12" it seems positive or neutral. For the other color, 
A(FA35W — F606W) ~ 0.0-0.2. Typical errors in these colors are 0.1 mags. 

A succinct description of the colors along the semi-minor axes is more difficult, as the 
errors are larger. Both colors are strongly affected by the scattering strip which produces 
a peak of emission at 3" in the NE, and 3" and 3." 7 at the SW. Outside these regions, 
A(FA35W — F814W) ~ 0.3-0.5 mags, with indications of a decrease beyond 5" along the 
SW side. 

The SBPs sample just narrow strips of all the data, and specially in regards to the 
colors may not give the full picture. Figures [T21 and [TBI show the A(FA35W — F606W) and 
A(FA35W — F81AW) colors of the disk relative to the star. We show the colors both before 
and after deconvolution. In both figures, the data have been convolved with a 4-pixel FWHM 
gaussian. Formally, the deconvolved image is the correct one to explore color information, 
but it is marred by streaks resulting from incorrect PSF matching, amplified by the Lucy- 
Richardson procedure. The effects of the scattering strip are manifested as two clumps both 
at the NE and SW sides of the star. As we mentioned before, the A(F435W - F814W) 
color is ~0.12 mags less beyond 2." 5 in the deconvolved data. 

These color maps make evident that the star is surrounded by an elongated nebulosity 
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8" in size (measured along the semi- major axes). Beyond this limit there is an abrupt change 
in the characteristics of the dust, as only a weak color signal is detected (however, nebulosity 
is detected beyond this limit). Within the mid-disk, A(F435W — F814W) increases from 
the <0.5 at the SW side, to >0.5 at the NE side, but it is mostly positive. If we assume that 
the SW side is closest to the Earth, this suggests a decrease in the value of the scattering 
factor g as a function of increasing wavelength. For A(FA35W — F606W) the situation is 
less clean-cut, although overall it is also true that the color index is larger along the NE 
than along the SW side. Comparing figures M and d2] indicates that the NE semi-minor axis 
is oriented along a particularly blue direction. 

In section 13.31 we noted that the arms have more neutral colors than the neighboring 
material. The effect is clearly seen in figure [131 on which the arms can be roughly traced by 
the darker color. Along the NW side of the disk, A(F435W — F&lAW) ~ 0.1 mags within 
3" , on the region occupied by structure 1. In the case of structures 2 and 3, the disk becomes 
neutral at the center of the features. (Along the smoothed SBPs the effect is only evident 
along the broad structure 2. See figure 0) The space within the arms has a non-neutral 
color in A(F435W - F8UW). 

4. Discussion 

As we have shown, these multicolor observations reveal a wealth of features in the cir- 
cumstellar environment of HD 100546. Any model of this system should be able to reproduce 
the following characteristics: 

• The circumstellar environment can be divided in three regions, with different brightness 
and color characteristics: the inner disk (from l."6 to 3"), the mid-disk (from 3" to 
8") and the extended envelope (beyond 8"). For all distances and at all angles, the 
disk is intrinsically brighter in the F814W band than in the other bands. 

• The inner disk is the region of the spiral arms. In addition to the well-known spiral 
structure, there are smaller arms never before reported. If we assume that the disk 
rotates counterclockwise the main arms trail the direction of disk rotation. Colors of 
the arms are more neutral than those of the surrounding material. The space between 
them becomes brighter as one goes to longer wavelengths, even after deconvolving by 
the coronagraphic PSF. The morphology does not change as a function of wavelength. 

• The mid-disk is brighter and redder to the SW than to the NE. If the larger brightness is 
interpreted as the result of scattering asymmetry, it implies that the SW side is oriented 
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towards the Earth. A description with a single value scattering asymmetry parameter 
implies g ~ 0.15, although there are local variations suggesting values as large as 
g ~ 0.2 The mid-disk has A(F435W - F8UW) ~ 0.5 and A(F435W - F606W) ~ 0.1 
although strong local variations are present. Comparison with NICMOS observations 
indicates that A(F435W - F160W) = 1.5 at 2." 5. The SW-NE color asymmetry 
suggests that g decreases with increasing wavelength. 

• An abrupt change of color and brightness marks the end of the mid-disk at 8" along the 
major axis. Farther out, we detect nebulosity to 12"-14". Overall, the shape is roughly 
consistent with what is expected from a flattened envelope with the same inclination 
as the inner disk. 

• Power-law fits to the SBPs produce indices steeper than —3 to 8". This dependence 
has also been observed in NICMOS and STIS data. 

A detailed model of this system is beyond the scope of this paper, but in the following 
sections we explore possible interpretations of these observations that may serve as a basis 
for such a model. We argue that the behavior of the SBPs with distance suggest that we 
are not detecting the optically thick disk responsible for the far-IR emission in the SED 
but only the optically thin envelope. The spiral arms are then structures in this envelope. 
Overall, the explanation for the colors of the dust remain a mystery, but we show that they 
are similar to colors of Kuiper Belt objects in the solar system. This suggests the processes 
like cosmic ray irradiation may be important to understand the dust color in HD 100546. 



4.1. The Color of the Circumstellar Material 

Given the large fractional infrared excess (Ldiskl L* = 0.51), we expect some part of the 
circumstellar material to be optically thick at visible wavelengths. This is not inconsistent 
with the very low relative scattered fluxes, as the amount of light scattered by an optically 
thick disk is controlled by the disk flaring angle. Neglecting the wavelength dependence of the 
scattering phase function, one expects that a situation in which every light path reaching the 
observer scatters off an optically thick structure to result in gray scattering. This i s because 



inters tellar dust albedo is approximately constant with wavelength in the optical (j Whitney 



19951 ). Therefore, the colors of the circumstellar material, especially those of material within 
«i3", which has been modeled as a flared, thick disk, are somewhat unexpected. On the 
other hand, optically thin lines of sight should result in blue scattering, if the dust around 
HD 100546 has the same optical constants and size distribution (wit h dust particl es ranging 



in size from tens of Angstroms to a few microns) of the diffuse ISM (lDrainell2003l ) 
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By ana lyzing the extinction as a function of distance for B-stars near the line of sight to 
HD 100546. iMalfait et al.l (119981 ) concluded that the reddening of the latter is consistent with 
the ISM column density at 100 pc. This means that the surrounding nebulosity does not 
contribute significantly to the reddening and the star is in front of the observed structures 
(or that the optical depth of circumstellar material in front of the star is ry < 0.05, from 
the error of the extinction measurements). In principle, some of the nebulosity could be 
associated with the nearby Lynds dark cloud DC 296.2-7.9, whose nominal center is only 
~25' away from the star. The ISM dust associated with the dark cloud should scatter more 
strongly in the blue than in the red. 

From figure [U the brightness in F814W along the major disk axis would have to be 60% 
smaller (or the brightness in F435W 60% larger) for the scattering to be gray (using 0.5 mags 
as F435W-F814W color). This is much larger than the uncertainties in the measured stellar 
colors (Section I2.3p . This red color is present even before the Lucy-Richardson deconvolution 
procedure, the effect of which is to reduce the F435W-F814W color at large angular distances. 
The disk is also red also in the F435W-F606W color, although not as strongly. 

Multicolor, resolved, scattered-light observations of disks around Herbig Ae/Be stars, 
that could be used to compare with our observations of HD 100546 are not very common. 
There are 14 Herbig Ae/Be stars for which the disk has been resolv ed^. Of those , six or seve n 



(in a ddition to HP 100546 these are: AB Aur, IGradv et al. 



2000; HD 150193A. lFukagawa et all2003l: PDS14 4N 



1999 



Perrin et al. 



HD 163296, IGradv et al. 



20061 : HD 142527. iFukagawa et al. 



20061 : and perhaps HD 169142. lHales et al.ll2006l ) have been resolved in scattered light. Be- 
sides HD 100546, HD 142527 (F6, 2 Myr-old) and PDS144 N (~A2, age uncertain) have 
observations in multiple bands - in the near-IR - performed by the same group, ensuring 
that the photometric calibration is consistent. In the case of HD 142527 those observations 
reveal a gray disk. PDS144 N is an edge-on system, and so the disk appears as a dark gray 
lane. Observati ons of T Tauri stars are m ore common, and for the most part reveal gray disks 
(e.g. TW Hya, IWeinberger et al.l 120021) . An interesting exception is GG Tau, which has a 
red, optically thick circum&man/ dis k (IKrist et al.l 2005b) . The large op tical color of the disk 
relative to the star (A(V — I c ) = 0.8. lKrist. Stapelfeldt. &: Watsonll2002l ) has been attributed 
to a combination of reddened illumination by the central starlight up on passing through in - 
ner circumstellar disks and the scattering characteristics of the dust (IDuchene et al.l 120041 ) . 
Observations of edge-on T Tauri stars reveal outflow-cavity walls blue in color (relative to 
their illuminatin g stars), and whose opt i cal characteristics can be exp lained by modified ISM 
extinction laws (jStapelfeldt et al.ll2003l : IWatson fc Stapelfeldtl 120041 ). 



3 Catalog of Resolved Circumstellar Disks, C. McCabe, www.circumstellardisks.org/ 
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If one considers the wavelength dependence of the phase fun c tion, it is possible to pro 



duce n on-grey, optically thick disks. iMcCabe. Duchene. fc Ghezl ( 120021 ) and iDuchene et al. 



(120041 ) argue that at least some fraction of the observed color of the GG Tau disk is due to a 
combination of the decrease of the scattering asymmetry pa rameter g with incre asing wave- 
length, and the inclination of the disk. If g is large (g ~ 0.5 in lMcCabe et al.ll2002l ). this effect 
can produce very red disk colors relative to the star, even in optically thick disks. Figures [121 
and [13] show evidence that the phase function plays a role in the color of the disk, because 
the dust becomes redder as the scattering angle becomes larger (assuming that the SW side 
is closest to Earth). However, if the angular dependence of the scattering function can be 
modeled by a Henyey-Greenstein function, an explanation based solely on the wavelength 
dependence of the phase function does not reproduce the colors. The Henyey-Greenstein 
function is such that if backward scattering results in A(FA35W — F81AW) = 0.5, forward 
scattering should result in A(FA35W — F814W) = —0.5. In other words, if the back of the 
disk is very red, the front should be very blue. This is not what is observed in HD 100546. 

Here we discuss three additional possibilities to explain the color of the observed circum- 
stellar material: (1) reddening due to dust very close to the star, (2) the dust characteristics 
of an optically thin envelope, and (3) selective absorption on the surface of the grains. 



4-1.1. Obscuring dust 

By analogy with GG Tau, we consider the possibility of obscuring material near the star 
which reddens the stellar flux illuminating the disk. Using the standard ISM extinction law 



with a normalized optical extinction value Ry =3.1 (ICardelli et al.lll989l ). we conclude that 
the 60% excess in F814W with respect to F435W, requires Ay = 0.65 mag to produce a gray 
disk. In other words, if there is a cloud of ISM material close to the star with Ay = 0.65 
mag, the disk would appear red although it would really be gray. Structures 1, 2, and 3 
would then be slightly blue. However, the SE Major axis F435W-F606W color is at most 
0.2 mags, which implies Ay = 0.45 mag. The same value of extinction can reproduce both 
colors if the extinction law changes more slowly in the optical. If Ry ~ 12, the colors will 
match with Ay m 1.3 mag. This would be in addition to the photometrically m easured 



stella r extinction (Ay=0.28 mag). Such large value of Ry has never been measured (IDraine 



2003|). 



Some stellar obscuration is expected from the models developed by iDullemond et al. 



( I200ll ). due to shadowing of the starlight by the inner disk rim. For the particular case of 
HD 100546, such a rim is quite small, which suggests that this source of optical depth is 
not important at large distances from what we call the inner disk. Notice, however, that we 
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detect red colors at all d i stance s >1."6, well within the ~4" disk radius predicted by the 
models by iDominik et al.l (120031 ) mentioned in the introduction. 



4-1.2. A Thin Envelope 



We may be detecting an optically thin envelope, as suggested by lVinkovic et al.l (120061 ). 
made up of a mixture of particles whose total opacity increases with increasing wavelength. 
This behavior of the opacity see ms to be common o f most debris disks with resolved circum- 
stellar material in scattered light (IKrist et al.ll2005al ). For debris disk systems the circumstel- 
lar material is more evolved (there is less gas, dust grains are larger, the dust is optically thin) 
than in Herbig Ae/Be stars. Observations of HD 141569, a 10 Myr old AO star surrounded 
by an optically thin debris disk s how that the disk co lors are A(F435W — F606W) = 0.1 



and A(FA35W - F8UW) = 0.25 (IClampin et al 



200 31) and recent observ ations of j3 Pictoris 



reveal similar colors along the spine of its disk (IGolimowski et al.l 120061 ) . No known debris 
disk has optical colors as large as HD 100546. 

A scattering opacity increasing to the red will result if the longest wavelength of ob- 
servation is comparable to the lower limit in the particle distribution, even for a standard 
astronomical silicate. The scattering efficiency is given by Q = a/na 2 , where a is the scatter- 
ing cross sectio n and a is the particle radius. Using the smoothed astrono mical silicate from 



Laor fc Drainei (119931 ). and integrating over the bandpass (see Eqn. 8 of IGolimowski et al. 



20061 ) we obtain the effective scattering efficiency, Q e s given the left panel of figure HU The 
secondary peaks have been smoothed, as they are the resu lt of the assumption o f perfectly 
spherical particles in the scattering opacities calculated by lLaor fc Drainei (119931 ) . 



Qeff has a broad maximum at a particle size comparable to the wavelength of observa- 
tion di vided by 4(n\ — 1), where n\ is the real part of the complex refraction index of the 
grain (lAugereau fc Papaloizoul 120041) . Noti ce that the Q e fr from figure HU is not the same 
quantity derived by IGolimowski et al.l (120061 ). In that work, the authors' "effective scattering 
efficiency" includes effects due to the dependence of the optical depth on the line of sight. 

To compare with observations, one should integrate the effective scattering cross section 
over the range of particle sizes, weighted by the particle size distribution. The resulting 
observable, the normalized surface brightness, is given by: 



I/F* oc 



,2-p 



Qefffg(0)da 



where / is the observed surface brightness, F* is the observed stellar flux and f g (9) 
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is the scattering phase function, which we assume to be a Henyey-Greenstein function. It 
depends on the particle size through g. We have assumed that the particle size distribution 
is dn/da oc a~ p . The proportionality constant depends on the surface density of particles as 
a function of distance from the star. 



The ratio of I/F* between two bands is shown in figure [14] (right), where we plot the 
A(F435Vr - F606W), A(F435W - F814W1 and the A(F435W - F160W) colors. For 



this figure we assume that p = 3.5 (jDohnanyi 



19691 ) and that the scattering angle is 9 =90 



degrees. In other words, these should be the colors along the major disk axis. 

If the lower limit in particle size is ~ 0.2 /im, colors like A(F435W — F606W) ~ 0.2 
mags and A(F435iy — F814W) « 0.5 mags are possible, but the predicted A(F435W — 
F1Q0W) color is too small, compared to the observations. If the lower limit in particle size 
is ~ 0.7 /im, A(FA35W - F160W) » 1.5 mags, A(FA35W - F606W) » 0.2 mags, but 
A(F435W — F814W) is too large. Analogo us results are obtained with pure graphite dust 
and with different grain porosities (as done in lGolimowski et al.ll2006h . We also experimented 
with different particle size distributions, ranging from uniform distributions (p = 0) to single- 
sized particles (p = oo). The three colors cannot be reproduced simultaneously. 

These arguments seek to explain the color of the scattered light as being a function of 
particle size, and assume that there is no selective absorption due to surface effects. The 
fact that the disk is not gray then implies the existence of "small" particles. For the dust 
model shown in the right panel of figure dH the disk becomes neutral in the optical if the 
minimum size particles are >5 /im. If the system does not have any gas, the lower particle 
size limit is set by radiation pressure. For HP 100546, a m %n is between ~5 and ~50 Atm , 



depending on the exact dust model used (ITakeuchi fc Artymowiczll200ll ; iLi fc Luninell2003l ). 
The presence of particles smaller than the radiatio n pressure limit is possible in systems 
with even modest amounts of gas (lArdila et al.ll2005T). In deed, any explanation of the richly- 
featured ISO spectrum near 10 /im (iMaffait et al 1998 ) require s the p resence of particles 
< 2/im in size close to the star. For example, Ivan Boekel et al.l (120051 ) models the 10 /im 
feature with dust whose crystalline component is dominated by forsterite grains 0.1 /im in 
radius. 

While models of the mid-IR spectrum consider dust much closer to the star than the 
one described in our observations, iBouwman et al.l (120031 ) argue that small forsterite grains 
(< 1/iin) are present in the disk at all distances. They attribute this result to collisional 
stirring by an inner (~ 10 AU) planet which expels small particles off the disk plane. The 
particles are then pushed out by radiation pressure tempered by gas drag. The lack of 
information about the gas content of the disk makes it impossible to predict the dynamics in 
detail. Segregation of dust by size is also possible if dust grains have grown and some have 
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settled in th e disk mid-plane. Smaller dus t particles will remain suspended at larger scale 



heights (e.g. iTanaka. Himeno. &: Idall2005l ) 



4-1.3. Selective Absorption 

If there is selective absorption (i.e. if there are strong albedo changes in the observed 
wavelength range, beyond what is expected from models based on ISM-like dust) the previous 
analysis will be incorrect. The scattering characteristics of the dust may then be dominated 
by the optical constants of materials deposited on the surfaces of the grains. In the solar 
system, Kuiper Belt and Centaur objects (KBOs) present relative col ors with respect to the 
sun ranging from neutral to A(B — I) 2.5 mags (ILuu fc Jewittlll996l ). This trend continues 
into the near-IR (IDelsanti et al.ll2004J ). Optical colors like those observed in HD 100546 are 
not uncommon (figure [T5l) . In this context, the F160W observation is within the range of 
colors sampled by the KBOs, having a normalized reflectivity (to V) of ~ 3. 



It has been argued (see lLuu fc Jewittlll996l ; IDelsanti et al.ll2004l and references therein) 
that the diversity of colors among KBOs and Centaurs is the result of the competing pro- 
cesses of formation of organic compounds on the icy surfaces (due to space weathering) and 
resurfacing (due to mutual collisions and cometary activity). The sources of space weather- 
ing are UV irradiation from the sun (most import ant in the case of solar flares), solar cosmic 
rays, and Galactic cosmic rays (iGil-Huttonl 120021 ) . Independently of whether or not this is 
the correct explanation for the observations presented here, weathering and resurfacing pro- 
cesses will also be at work around HD 100546, because the water-ice sublimation boundary 
(T~120 K) occurs at R ~ 30 AU for large particles. 



Delsanti et al.l (120041 ) developed a model of the re ddening process in KBOs, in which 



the reddening timescale (r ~ 10 yrs. IShul'Man 



19721) is controlled by the Galactic cosmic 



ray flux. In the case of j3 Pic, iGolimowski et al.l (120061 ) used this model (with the timescale 
scaled by the stellar luminosity) to conclude that the dust grains escape the system before 
they are significantly altered. However, this scaled model may not be applicable to (3 Pic or 
HD 100546. As they are early type stars, their astrosphere will be smaller than the sun's 
(reducing the role of stellar cosmic rays) and the effect of galactic cosmic rays will be more 
important. Furthermore, for HD 100546 the UV flux (A < 1500 A) at 300 AU is ~300 
times larger than the solar UV flux at 40 AU. Recent analysis of the colors of the debris 
disk around HR 4796 (an AO star, « lOMyrs old) also suggest a scattering opacity similar 
to that of some KBOs (A. Weinberger, 2006, personal communication). In summary, space 
weathering is likely to play a role in our understanding of the colors of dust around A stars. 
However, a model appropriate for the conditions around those stars has not been developed. 
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4- 1-4- Summary 

Our observations reveal two separate, independent facts. First, the circumstellar mate- 
rial has a non-zero color relative to the star. This is true even within 3" , where the thermal 
emission has been modeled as being due to an optically thick disk. The color suggests that 
there is scattering asymmetry in the dust, extinction close to the star, optically thin lines 
of sight, or selective absorption. Second, the circumstellar material is red, even to large dis- 
tances of the star, suggesting that all has been processed in a similar way. As we have shown, 
extinction requires an anomalously large value of Ry to explain both sets of colors: values of 
Ry this large have not been observed. An optically thin disk halo composed of astronomical 
silicate cannot explain the measured scattering asymmetry. The colors are within the range 
of colors observed in KBOs, which are believed to be produced by selective absorption on 
the surfaces of dust gains. 

As in GG Tau, it is possible that more than one process is at work in the system. 
Beyond 3", the role on inner disk extinction on the light scattered off the mid-disk should 
be relevant for our understanding of the system. The intense UV radiation coupled with the 
effect of Galactic cosmic rays may redden the dust before is dissipated. These ideas await 
better modeling efforts. 

Except for an explanation of the colors based only on extinction, the remaining ex- 
planations imply that the dust opacity increases with wavelength. The explanations based 
on optically thin scattering of a yet-to-be-determined dust mixture and that of a KBO-like 
material are not necessarily different, as the polymerized ices believed to be responsible for 
the KBO colors may provide the required dust characteristics. Formally, dust coated with 
this kind of ice mixture will be red, no matter what its optical depth. However, as we show 
below, there are other pieces of evidence that suggest that we are observing an optically thin 
dust distribution. 



4.2. The Envelope 



Our fit to the SBPs within 2." 7 using the lGrady et al.l (120011 ) ranges (see bottom portion 
of table [5]) is consistent with a r~ 3 dependence, which suggests that t he structure responsible 
for th e scattered light is a geometrically thin, opticall y thick disk ([Whitney fc Hartmann 
19921). This same r~ 3 depe ndence has been observed by iGrady et al.l (120011 ) (from 1 to 2.7") 
and lAugereau et al.l (120011 ) (from 0.5 to 2.5"). However, the SBP for a flat disk is only 
10~ 7 / (ft 2 arcsec -2 , where is the field distance (in arcseconds) to the star (jWhitney fc Hartmann 



19921 ). This is about 3 orders of magnitude less than what is observed. A flared, optically 
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taeawa et al. 


2004) 


Gradv et al. 


(1999) 



light surface brightness profiles decrease with stellar distance as r~ 3 
This near-IR dependence is steeper than the r -2 dependence found by 
in the optical. 

On the other hand, the SBP of an optically thin disk may present any radial power-law, 
depending on the surface density of the scatterers and/or changes in the scattering cross 
section as a function of distance. If the material is optically thin along the line of sight, the 
surface brigh tness profile then im plies a surface density of r _1 or r~ L8 (depending on whether 
one uses the iGrady et al.l (120011 ) tracks or our 0."25 wide tracks), ignoring variations in the 
scattering cross section. However, it is difficult to understand how the material around 
HD 100546 can be optically thin within 3" , because the dust in this region is the bulk source 
of the IR emission. A similar power law dependence is observed from 0."5 to 3" suggesting 
that the disk does not change from optically thick to optically thin within this range. It is 
unlikely that the transition to optically thin occurs within 0."5 (50 AU) from the star as in 
this case, the spectra l signature from the d isk interior would appear much hotter than the 



50K observed (e.g. iBouwman et al.ll2003l ) 



An argument analogous to this led IVinkovic et al.l (120061 ) to posit the existence of a 
flattened optically thin envelope or halo surrounding a geometrically thin, optically thick 
disk. The disk surface (the surface where the optical depth from the observer to the disk 
« 1) has to be shallow enough that it does not intercept sig nificant amount s of th e stellar 
radiation, but otherwise is unconstrained. The objective of IVinkovic et al.l (120061 ) was to 
explain the 3 /im bump mentioned in the introduction. They concluded that the envelope 
has constant optical depth in the visible given by ry = 0.35, although this refers only 
to relatively hot material responsible for the bump. In the context of this picture, our 
observations would only be detecting light scattered by the optically thin component. The 
optically thick disk is undet ected in thes e optic al observations. A remarkably similar picture 
has also been suggested by lMeeus et al.l (120011 ) (see their figure 8, top panel) to explain the 
mid-IR SED. 

In summary, the distance dependence of the SBPs suggest the presence of an optically 
thick shallow disk with a geometrica lly th ick, optically thin, fla ttened envelope, similar to 
that advanced by iMeeus et al.l (120011 ) and IVinkovic et al.l (120061 ). No thermal or scattering 
model with these characteristics has been published and so it is not clear if this picture may 
reproduce the SED and the lar ge fractional infrared excess. A flat, optically thic k disk can 
only produce L IR /L* = 0.25 (I Adams fc Shul 1 19861 ; I Whitney fc Hartmannl Il992l ). A more 
realistic, flared, optically thick disk, would need to have a total opening angle of 60° to 
explain the excess. If the disk is flatter than this, the envelope would have to cover a larger 
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angle (as seen from the star) and/or have a large radial optical depth. Our observations 
show an envelope extending to 1200 AU from the star, but whether this will be enough to 
reproduce the SED remains uncertain. 



4.3. Spiral Arms or Shadows? 



Various explanations have been advanced to understand the spiral structure seen in the 
scattered-light images. On the one hand, it is possible that they are "true" arms, due to 
spiral density waves or tidal effects. iQuillen et al.l (120051 ) conclude that the disk would have 
to be unreasonably thin for the Toomre Q parameter to be less than one, suggesting that 
the disk is stable to spiral-density waves. 



Guillen et al.l (120051 ) also examined the possibility that the spiral structures are driven 



by a compact object external to the disk, either a bound planet or a passing star. They per- 
formed hydrodynamical simulations (i.e. a pure gas disk) of the effect of bound planets and 
unbound stellar encounters on a 300 AU disk and concluded that a 20 Mj mass planet/brown 
dwarf in an eccentric orbit would be necessary to excite the observed spiral structure, which 
they argued would be detectable in the NICMOS images. They also concluded that the 
hypothesis of an encounter with a passing star lacks suitable candidates. 



Could we detect such a body in our images? Figure [161 shows the 3a detection limits for 
point sources around HD 100546. Because the differences in detection along the major and 
minor axes are small, we average them. To obtain these curves the number of c ounts has been 



transf ormed into instrumental vegamag magnitudes using the quantities from ISirianni et al. 



(120051 ). Such transformation ignores c olor terms. Let us c onsider a 10 Mj mass planet, 10 
Myrs old. According to the models by iBaraffe et al. such an object has Mj c = 14.8, 

which implies mj c = 19.8 at 100 pc, and mp8i4w ~ 19.4, in the instrumenta l vegamag 
system. The V-I color of such an object would be 5.5 mags (IBaraffe et al.l 120031 ) . and so it 
would be undetected in the F606W band. We searched for objects that were present in the 
F814W images but not in the F606W ones, but did not find any. Clearly, we would have been 
able to detect such object in a ny place outside the mask. This planet would have m# w 16.3 
mags. lAugereau et al.l (120011 ) do not provide detection limits as a function of distance for 
their observations, and s o it is difficult to say if it would be visible in the NICMOS images. 
The objects predicted by lQuillen et al.l (120051 ) to generate the structures would certainly be 
detectable in the images presented here. 

A number of researchers (e.g. iBouwman et al. I l2003l : iGrady et al. 1 120051 : lAcke fc van den Anckerl 

20061 ) have indicated that the presence of a hole in the gas and dust distribution at w 10 AU 
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suggests the pres ence of a planet. In this case, the perturber is internal to the disk, unlike 
the simulation by lQuillen et al.l (120051 ). The effects of this configuration on HD 100546 have 
not been explored. 



Warps in the cast shadows on the outer disk that may appear as spiral arms. iQuillen 



(120061 ) modeled the effect of warps on a disk that is optically thin perpendicular to the disk 
plane (within 80 AU) but partially optically thick along it. In this model, the dark lanes 
are due to opacity of the mid-plane and the bright lanes represent optically thin material 
partially oriented along the line of sight. The model predicts that, as the disk is observed 
at wavelengths in which the opacity is larger, the dark lanes will change little, or they may 
become darker (because the mid-plane optical depth will increase) and the bright lanes will 
be brighter (as the optical depth along the line of sight increases). 

If we assume that the dust opacity increases from F4 3 5W t o F814W, the brightening 
of the inter-arm space contradicts the models from iQuillenl (120061 ) and is more in tune with 
what one expects if the structures are "real arms": changes in the volumetric density of 
material separated from the inner disk by relatively empty regions, that become brighter at 
optically thicker wavelengths. However, as mentioned in section I4.2[ we believe the spiral 
arms are observed in the envelope, not in the optically thick disk. They may well reflect 
structures in the optically thick disk, but the latter is undetected. 

In summary, the observations presented here are inconsistent with the published simu- 
lations that seek to explain the origin of the spiral arms. Whether a model can be developed 
that generates the observed structure with an interior planet, with a very small exterior 
planet, or via warps in a more realistic, larger disk, remains uncertain. 



5. Conclusions 

We present ACS/HRC coronagraphic images of HD 100546's circumstellar environment 
in the F435W, F606W, and F814W bands. The star is a B9.5, >10 Myr old, 103 pc away 
from the sun. The observations were performed with the 0."9 occulting spot. For each band, 
we observed the system in two different telescope rolls. To improve the contrast between the 
star and the circumstellar material, we subtracted a PSF reference star from HD 100546. 
We used HD 129433, which has a similar color as the target. To scale the two, we used direct 
measurements in the F435W band, coupled with color measurements of HD 100546 taken 
from the literature. We deconvolved the images using a simulated stellar PSF, obtained with 
the PSF simulator Tiny Tim. Both the deconvolved and non-deconvolved images in each 
band are available from the Journal Web page. 
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Subtraction residuals dominate the dataset within ~1."6. Color-combined images re- 
veal a large scale structure of scattered light to distances «14" away from the star. The 
well-known "spiral arms" (which we call structures 1, 2, and 3) of the disk are clearly seen. 
Photometric analyses by other groups reveal that the star is in front to the observed nebu- 
losity. 

The scattered-light relative fluxes measured from l."6 to 10" in the F435W, F606W, 
and F814W bands are (L circum / L*) A = 1.3±0.2 x 1(T 3 , 1.5±0.2 x 1(T 3 , and2.2±0.2 xl(T 3 , 
respectively. If we interpret the circumstellar nebulosity within 2" as due to an inclined disk, 
the PA of the semi-major axis is ~ 145 degrees, with an inclination of ~ 42 degrees. 

We analyzed the scattered light by taking median cuts along segments 0."25 wide, 
oriented along the disk's principal axes. For these profiles, the 3a detection limits are 
2 x 10~ 7 /arcsec 2 (corresponding to Vv eg a=23.4 magnitudes/arcsec 2 ). Analysis of these 
profiles reveal that the circumstellar material scatters more strongly towards larger wave- 
lengths. Along the major axis, the A(F435W — F606W) color ranges from to ~ 0.2. The 
A(F435W — F814W) varies between m 0.5 and ~1, increasing slowly with distance. Com- 
parison with NICMOS results show that A(FA35W - F1Q0W) ml. 5. The SW side of the 
nebulosity is brighter (at all distances) than the NE side, by factors between 1 and 3, which 
we interpret as being due to forward scattering. This suggests that the SW side is the side 
closest to the Earth. The brightness ratio can be explained invoking a Henyey-Greenstein 
scattering phase function with g = 0.15 (assuming an optically thin disk). However, there 
are local variations that suggests values of g as large as 0.23. Evidence of a change in g with 
wavelength is deduced by the fact the the SW side of the disk is bluer than the NE side. 

Conceptually, we divide the disk in three regions, described by their size as measured 
along the major axis. The inner disk (from l."6 to 3" ) has a size comparab le to the optically 
thick disk used to explain the system's far-IR SED ( Dominik et al.lfeood ). This is also the 



region occupied by the spiral arms. The outer edge of the mid-disk (from 3" to 8") is 
marked by a factor of two decrease in brightness and by the disk color becoming almost 
neutral. Beyond that we detect nebulosity to ?^14" along the major axis. The SBPs are 
consistent with the mid-disk and the extended nebulosity having the same inclination as the 
inner disk. 

The well-known spiral arms are clearly detected in the images. They are seen as a 
narrow bands of nebulosity that wrap around the disk. If the SW side is closest to Earth, 
spectroscopic observations that spiral structure trails the rotation. There are no significant 
morphological differences among the images in the different bands. The space between the 
arms and the part of the inner disk closer to the star becomes brighter at longer wavelengths, 
but the peak brightness of the structures does not change much, i.e. the structures become 
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grayer at longer wavelengths. The brightening of the inter-arm space with wavelengths is 
inconsistent with models that assume that they are due to the effects of a warped disk. We 
deproject the disk and subtract a smoothed version of itself to increase the contrast of the 
sharp structures. This procedure suggest the existence of smaller, fainter arms, not easily 
seen in the original images. The NW arm along the major axis is split in two and the SW 
arm is revealed to have a companion along the NE side. The strongest arms seem to trace 
a four- armed spiral. 

Dynamical models that seek to explain the structures as due to perturbations from 
external planets predict the presence of objects that are not seen in the images presented 
here. Observations of hot gas suggest the presence of a planet at 10 AU, but detailed models 
of the effect of this object in the disk of HD 100546 have not been published. The cause of 
these structures and the mechanism for their preservation remains unknown. 

If one ignores the wavelength dependence of the scattering phase function, the fact that 
the inner disk has a color is surprising, given that the optical albedo derived for ISM grains 
is constant with wavelength, and that a large optical depth is predicted for this part of the 
disk. Most young stellar objects in a similar evolutionary stage are gray. The more evolved 
debris disks have mostly red colors, although not as red as those of HD 100546. If the colors 
were due purely to the behavior of the scattering phase function with wavelength, the disk 
would be bluer in the front side that what is observed. We examine the possibility that the 
color is due to the presence of dust between the star and the disk, beneath the coronagraphic 
mask. This would require different amounts of dust to explain different colors among the 
same direction, if the dust is ISM-like. Ry would have to be « 12 before we can reproduce 
the colors with the same amount of dust. We also explore the effect of the minimum particle 
size on the colors of a mixture of "astronomical silicate" grains. We conclude that a size 
distribution dn/da oc a~ p with p = 3.5 is not capable of reproducing the colors. This general 
conclusion is also true for pure graphite grains, non-compact silicate grains, and for all values 
of p. In the context of an explanation for the color based on particle size, the color of the 
system implies the presence of particles <5/im at distances up to ~800 AU from the star. 

We show that the observed colors are within the range of colors observed in KBOs in 
the solar system. For KBOs the colors are believed to be the result of a competition between 
space weathering (cosmic ray and UV flux irradiation) and resurfacing due to collisions. This 
is a promising scenario and deserves to be explored more, by observing the circumstellar 
material at other near-IR bands. 

Within 2". 5, the major axis SBPs has a steep r~ 3 ' 8 dependence. This suggests that the 
structure responsible for the material is not an optically thick disk. A flat optically thick 
disk would result in SBPs with a r~ 3 dependence on distance. A flared disk would result in a 
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shallower power law. The distance dependence of SBPs requires invoking the presence of an 
optically thin halo or envelope in addition to the o ptically thick disk. Our o bservations are 



then c onsistent with the qualitative pictures from iMeeus et al.l (120011 ) and IVinkovic et al. 



( 120061 ). in which an optically thick disk is responsible for the far-IR emission but a large 
scale- height, optically thin disk (with a scale-height larger than predicted by standard disk 
models) surrounds the optically thick component. Given its colors, the envelope cannot 
be composed of remnant infalling ISM material but has to have been reprocessed by the 
star+disk system. The optically thick disk itself is undetected, suggesting that it has a very 
shallow surface profile. More detailed models are required before the characteristics of this 
envelope are better known. 

The crowded field of HD 100546 suggests a way to probe the envelope material, by taking 
spectra of stars 3, 4, 5, 8, and 9 (Figure ED- However, for optically thin dust, the extinction 
in V due to the envelope is of the order of SB x 47T0 2 , where SB is the normalized surface 
brightness at F606W, and <fi is the angular distance between HD 100546 and the measurement 
position. This implies that at 5", Ay ~ 0.001 mags. While this is a lower limit (some light 
may be blocked by the inner disk) trying to detect a small amount of extinction on the 
spectrum of a star 5" away from a source forty thousand times brighter (comparing star 4 
to HD 100546), is a formidable obstacle for this kind of experiment. 

Further observations and models are required before we can fully understand this sys- 
tem. In particular, it is crucial to measure the cool gas content of the circumstellar material, 



which will determine the dynamics of the dust (lArdila et al.ll2005l ). Additional coronagraphic 



observations at other NICMOS bands will provide much needed color points at longer wave- 
lengths, which will help decide whether or not a new kind of dust model is needed to un- 
derstand old protoplanetary (and/or young debris) disks. The observations presented here 
show how powerful multicolor coronagraphic data can be in advancing the understanding of 
protoplanetary disk evolution. 
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Table 1. Log of Coronagraphic Exposures 



Object 


Band 


Obs. Start 


Total Exp.Time 


P.A. a 


Gain 






(UT) 


(sees) 


(deg.) 


(e~ /photon) 



HD 100546 



HD 129433 b 



F435W 


2004-04-26 13:51:34 


3 


-138.31 


4 




2004-04-26 22:39:00 


160 


-138.35 


2 




2004-04-26 22:43:00 


2400 


-138.35 


2 




2004-04-26 23:54:00 


160 


-118.35 


2 




2004-04-26 23:58:00 


2530 


-118.35 


2 


F606W 


2003-03-26 19:42:39 


3 


179.09 


2 




2003-03-26 03:53:00 


130 


179.05 


2 




2003-03-26 03:57:00 


2600 


179.05 


2 




2003-03-26 07:05:00 


130 


-152.95 


2 




2003-03-26 07:09:00 


2600 


-152.95 


2 


F814W 


2003-03-26 19:42:39 


6 


179.09 


2 




2003-03-26 02:25:00 


130 


179.05 


2 




2003-03-26 02:28:00 


2350 


179.05 


2 




2003-03-26 05:33:00 


160 


-152.95 


2 




2003-03-26 05:37:00 


2520 


-152.95 


2 


F435W 


2004-04-27 13:56:20 


2 


151.52 


4 




2004-04-27 01:32:00 


16 


151.48 


2 




2004-04-27 01:34:00 


160 


151.48 


2 




2004-04-27 01:40:00 


400 


151.48 


2 




2004-04-27 01:48:00 


1300 


151.48 


2 


F606W 


2004-04-13 19:45:02 


1 


117.76 


2 




2003-03-26 09:33:00 


16 


117.72 


2 




2003-03-26 09:34:00 


1020 


117.72 


2 


F814W 


2003-03-26 19:45:11 


2 


117.72 


2 




2003-03-26 09:13:00 


16 


117.72 


2 



a Position angle of image y axis (deg. E of N) 
b PSF Reference star for HD 100546 



Note. - 



Only observations made with the parameter GAIN=4 preserve the number of 



counts in a saturated imag 
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Table 2. Photometry 



HD 100546 (B9.5Vne) HD 129433 (B9.5V) 



F435W 



e-/sec (10 7 ) 2.5 (3%) 6.4 (3%) 

Flux Density (Jy) 8.3 (3%) 21.2 (3%) 

Bvega a 6.70 (0.03) 5.69 (0.03) 



F606W 



e"/sec (10 7 ) 4.9 (4%) 11.1 (3%) 

Flux Density (Jy) 7.3 (4%) 16.4 (3%) 

Vvcga a 6.69 (0.04) 5.79 (0.03) 



F814W 



e'/sec (10 7 ) 1.9 (5%) 4.1 (3%) 

Flux Density (Jy) 5.4 (5%) 11.5 (3%) 

Ivega a 6.64 (0.05) 5.82 (0.03) 



a Vegamag system. The instrumental system produces the same 
magnitudes within the errors. 



Note. — HD 129433 is the star used a PSF reference. For F435W 
the count rate (e~/sec) is measured from the image. The F606W 
and F814W count rates are estimated from the F435W values and 
synthetic photometry using proxy spectra. Th e F435W measure- 
ments have errors of ~3% ( ISirianni et al.ll2005l ). For HD 100546, 
the errors in F606W and F814W include errors in the measured col- 
ors of the target star in these bands, while the errors in HD 129433 
are only those of the measured photometry in F435W. 



Table 3. Standard Photometry of Field Stars 



Number 


B 


V 


I 


1 


23.7 ±0.5 


22.2 ±0.1 


20.8 ±0.1 


2 


> 24.5 


24.2 ±0.3 


22.1 ±0.3 


3 


22.1 ±0.2 


21.2 ±0.09 


20.2 ±0.1 


4 


19.00 ±0.05 


18.11 ±0.05 


17.14 ±0.05 


5 


19.37 ±0.06 


18.42 ±0.03 


17.44 ±0.05 


6 


22.1 ±0.2 


21.8 ±0.1 


21.1 ±0.2 


7 


> 24.0 


23.8 ±0.3 


22.5 ±0.3 


8 


23.2 ±0.4 


21.3 ±0.1 


19.00 ±0.06 


9 


19.5 ±0.06 


18.60 ±0.05 


17.64 ±0.05 



Note. - Aperture photometry for field stars, ob- 
tained from non- deconvolved images. See figure [5] for 
a key to the numbers. Color transformations from 
the instrumen t al to the standard system given by 



Sirianni et al.l (120051 ) . The error includes uncertain- 



ties in the measurement of the local sky and Poisson 
noise. It does not include the (systematic) errors intro- 
duced by transformation to the standard photometric 
system. 
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Table 4. Geometric Parameters for the Disk 



Reference 


Inclination (cleg from face-on) 


PA of Major Axis 


This Work 


42±5 a 


145 ±5 


Pantin et al. (2000) 


50 ±5 


127 ±5 


Grady et al.(2001) 


49 ±4 


127 ±5 


Augereau et al. (2001) 


51 ±3 


161 ±5 



a From isophot fitting between l."6 and 2" from the star. The numbers 
quoted are averages for all bands. This determination ignores the angular 
dependence of dust scattering. 
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Table 5. Power-law fits to the SE semi- major axis 





SE Maj 


or Axis 




Distance 


F435W 


F606W 


F814W 


l."6-5."0 


-3.8 ±0.1 


-3.8 ±0.1 


-3.7 ±0.1 


5."0-8."0 


-3.4 ±0.3 


-2.8 ±0.3 


-2.7 ±0.3 


8."0-12."0 


-0.6 ±0.4 


-0.4 ±0.3 


-1.2 ±0.5 




Gradv ct al. ( 


20011 ranges 




l."6-2."7 


-3.6 ±0.4 


-3.5 ±0.4 


-3.5 ±0.3 


2."7-5."0 


-4.2 ±0.2 


-4.2 ±0.2 


-3.9 ±0.2 


5."0-8."0 


-2.5 ±0.3 


-2.5 ±0.3 


-3.0 ±0.3 



Note. — The "SE Major Axis" fits are obtained 
along a strip 0." 25- wide, centered on the s tar and 
with PA=145 deg. For the lGradv et al.l l|200lt ) ranges, 
PA=127 deg, each strip is 0."46 wide and the fits are 
performed on the average of the SE and NW axes 
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Fig. 1. — Surface brightness maps of the circumstellar environment of HD 100546, in a loga- 
rithmic stretch. All images have been normalized to the stellar brightness in their respective 
bands. The approximate size of the coronagraphic mask is shown as a black circle 1" in ra- 
dius. The top row has a different color stretch and spatial scale than the other two, in order 
to showcase different regions. The top stretch goes from 10 -7 arcsec -2 to 10 -3 arcsec -2 . The 
middle and bottom rows show a stretch from 10 -5 arcsec -2 to 10 -3 arcsec -2 . The contours 
in the bottom row are obtained from images that have been heavily smoothed. The contour 
values are 2, 3, 4, 5, 6, 7, 9, 12, and 15 xlO -5 arcsec -2 . 
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SE Semi — Major axis 




1 10 
Distance from star (arcsec) 



Fig. 2. — The effect of the deconvolution in the surface brightness profiles. The plot shows 
the median surface brightness in segments 0."25 long, taken along the SE extension of the 
major axis of the disk (145 degrees East of North). Each profile has been smoothed by a 
median boxcar 10 pixels wide, and offset by an arbitrary amount in the y-axis. The dashed 
lines show the profiles before, and the solid lines show the profiles after the deconvolution. 
The dotted line is placed at l."6 from the star, the distance within which the images are 
dominated by subtraction residuals. Note that the effect of the deconvolution is to sharpen 
the features within ~2."5 At larger distances, the deconvolution is most important for the 
F814W band, resulting in a surface brightness profile reduced by ~15% with respect to the 
non-deconvolved one. 
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Fig. 3. — PSF-deconvolved images of the circumstellar environment of HD 100546. The 
color and spatial scales are the same as shown in Figure [TJ In order to draw the contours, 
the images have been smoothed like those in the bottom row of figure 1, and then convolved 
with a 10-pixel-FWHM gaussian (see Section 2.6). 



Fig. 4. — RGB composite of the individual images. As in previous figures, the black central 
circle is 1" in radius. North is up and East is left. Before combining them, each image has 
been normalized by the stellar brightness. In order to map the number of counts to the 
RGB channels, we have used the algorithm described by iLupton et al.l (12004 ). Left: Color 
combined image. Right: Color combination of the deconvolved images. 
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Fig. 5.- 



Guide to table [3] (Color-combined image of non-deconvolved bands). 
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Fig. 6. — The inner circumstellar region as seen in the F435W band. We have annotated 
structures 1, 2, and 3, which are observed in all bands. The scaling is the same as in the 
mid-panels of figure [TJ If the SW side is closes to Earth, spectroscopic observations suggest 
that the disk is rotating counterclockwise. 
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Fig. 7. — Unsharp masking of the non-deconvolved image in the F435W band. Upper left: 
Deprojected image, with each pixel corrected by scattering asymmetry (g = 0.15), and 
multiplied by <ft 2 , where <fi is the angular distance to the star. Notice that the hemispheric 
subtraction residuals are circular in the original image, so they look elliptical here. Upper 
right: The same image, smoothed by a gaussian kernel with 30 pixels FWHM. Bottom 
left: Unsharp masking, the result of subtracting the upper-right image from the upper- left. 
Bottom right: Same as bottom left, with features identified. We only mark those features 
that are identified in all the bands. Features la, 2a, and 3a correspond to structures 1, 2, 
and 3 in figure El 
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Fig. 8. — Top: Surface brightness profiles (SBPs) of the deconvolved images divided by the 
stellar flux. Bottom: Disk colors relative to the star. The SBPs are calculated by taking the 
median flux within segments 10 by 1 pixels (0."25 by 0."025) along an angle 145° East of 
North (the PA of the disk), to 13" away from the star. In the plot, the profiles have been 
also smoothed using a median boxcar 10 pixels long. The dotted lines mark l."6, 5", 8", 
and 12", corresponding to the circles shown in figure [10] and the ranges for the power- law fits 
(table EJ). The errors are shown only at representative points, as filled circles. On the scale of 
the top panels they are of the order of the size of the circles. Subtraction residuals dominate 
within l."6. In the top right panel, the squares are STIS measurem ents (jGradv et al.l 12 001) 
and the triangles are NICMOS measurements in the F160W band (jAugereau et al.ll200l[ ). 
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Fig. 9. — Same as in figure [8] but centered at an angle 55° East of North. Between 2." 5 and 
3." 7 in the NE direction, we see the effect of the scattering strip as an increase in the colors. 
Beyond 11." 75 in the NE, we reach the end of the region sampled by both rolls in F435W. 
Beyond 8" in the SW the profiles are affected by the large coronagraphic spots. 
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Fig. 10. — Annotated, deconvolved image of F435W divided by the stellar flux, to use as a 
guide interpreting the surface brightness profiles of figures |8] and [9j The circles indicate radii 
l."6, 5", 8", and 12", corresponding to the limits of the power-law fits. The rectangular 
boxes are 0."25 wide, and are oriented in the position angles extracted to produce the SBPs. 
The image is in a logarithmic stretch, with limits 5 x 10 -6 arcsec -2 to 10 -3 arcsec -2 . 
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Fig. 11. — Comparison of the profiles at both sides of the disk. The profiles of different 
bands have been shifted by an arbitrary amount. Blue, green, and red colors correspond to 
F435W, F606W, and F814W respectively. (Left) Solid: SE side; Dashed: NW side. (Right) 
Solid: NE side; Dashed: SW side. Also shown are the power-law fits for the F435W band in 
the major axis (starting closest to the star, the indices are -3.8, -3.4, and -0.6) 



Fig. 12.— The A(F435W - F606W) = -2.5 \og(SBP F435W / SBP F606W ) color. Left: Color 
image obtained from non-deconvolved data. Right: Color image obtained from deconvolved 
data. In the NE the effect of the scattering strip is seen as a brightening of the colors. These 
features have weaker counterparts on the SW. 




Fig. 13.— Same as figure [13] 

-2.5 \og(SBP F435W /SBP F8UW ). 



but with A(F435W 



- F8UW) 
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Fig. 14. — The color of the dust. Left: The scattering efficiency (averaged over observation 
band) as a function of grain size for the standard astronomical silicate. Blue (solid): F435W; 
green (dotted): F606W; red (dashed): F814W; black (dot-dash): F160W. Right: Expected 
color of the dust as a function of minimum particle size. We assume a size distribution in 
the form of a" p , where p=3.5. Blue: A(F435W - F606W). Red: A(F435W - F814W). 
Black: A(F435W - F160W). 
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Fig. 15. — Colors of KBO and Centaurs (black lines), compared with the measured colors of 
HD 100546. The ordi nate axis is the norm alized reflectivity, Rn = 10°' 4A ( mA_ ^, where V is 
the star's magnitude (ILuu fc Jewittlll996l ). For the HD 100546 we use the colors measured 
3.5" away from the star along the SE semi-major axis, A(F606W — FA35W) = —0.18, 
A(F606W-F8UW) = 0.29, and A(F606W -F160W) = 1.24. The colors were transformed 
from the instrumental to the standard vegamag system. The data for the solar system objects 
was taken from table 2 of iDelsanti et al.l (120041 ) . We include only those objects with BVRIJH 
photometry. Notice that the colors of the circumstellar material are well within the range of 
KBO colors. 
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Fig. 16. — Detection limits for point sources. The curves show the 3a detection limits for 
point sources in F606W (dashed) and F814W (solid) magnitudes. The curves were obtained 
by planting stellar sources at different distances from the star and calculating how bright 
would those sources have to be to be detected at the 3a level. The noise is calculated based 
on the local background value. These detection limits are averages of detection limits along 
the same NE, SE, SW, and NW directions used to obtained the profiles from figures M and 
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